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IMPROVED ULTRASONIC VOLUMETRIC IMAGING BY 
COORDINATION OF ACOUSTIC SAMPLING RESOLUTION, 
VOLUMETRIC LINE DENSFTY, AND VOLUME IMAGING RATE 

5 This invention relates to ultrasonic diagnostic imaging and, more 

particularly, to controlling the relationship of the acoustic sampling resolution, the 
desired output line density, and the volume imaging rate in ultrasonic volumetric 
imaging systems. 

Ultrasonic diagnostic imaging systems are now capable of scanning 

10 volumetric regions of a body for the production of three dimensional images of the 
volumetric region. Since many more beams are necessary to scan a volumetric region 
as compared to the planar region of a two dunensional image, the time required to scan 
a voktmetric region can be much greater, causing the rate at which volumetric images 
are created to be relatively low. One approach to maintaining an acceptable image rate 

15 is to predetermine a constant number of transmit beams which will be used to scan a 
nominal volumetric region for a given procedure such as cardiac imaging. As the user 
adjusts die depth of the image field to encon^ass a greater depth than that of the 
nominal volume, the frame rate will decrease as greater time is required to receive 
echoes from the greater depth. If the user adjusts a lateral dimension of the nominal 

20 volume so that a wider volumetric region is scanned, the transmit beams are spread out 
more widely to scan the wider volume and the beam density declines. This decline in 
beam density can result in a spatial undersampling of the volumetric region as the beam 
density decreases. For some applications a minimal spatial undersampling of &e image 
volume may be hardly noticeable. However, for other applications deleterious image 

25 artifacts will appear. Spatial undersampling of a planar or volumetric region will give 
rise to a shinmiering effect in the image, and it may seem as tihtough the image is being 
viewed through a grate or screen. In certain diagnostic applications such as a search for 
lesions of the liver, the pathology is often diagnosed by discerning subtle variations of 
the texture of the liver in the image. The speckle pattern of the ultrasound image can 

30 play a role in this diagnosis as the clinician looks for subtle change in the speckle 
pattern of the imhgp of the liver. Such subfle differences can be masked by tfie 
scintillating or shimmering arti&cts of spatial undersampling. Accordingly it is 
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desirable to prevent or at least control spatial sampling arti&cts so that such diagnosis 

will not be impeded. 

In accordance with the principles of the present invention, an ultrasonic 

volumetric imaging system is described in which spatial sampling is controlled by 
5 control of the acoustic imaging point spread function. In an illustrated embodiment the 

acoustic imaging point spread function is coordinated witti tiie line density of &e 

volumetric region to produce a desired spatial sampling of the volumetric region. 

Through such control an accq)tabl6 level of spatial sampling artifacts may be 

maintained as the size or shape of flie volumetric region is changed. In accordance with 
10 another embodunent of the invention, the scanning of greater depths may be afforded by 

control of the point spread fiinction within acceptable levels of acoustic output 
In the drawings: 

FIGURE 1 illustrates an idealized beam intensity in one dimension. 

FIGURE 2 illustrates the idealized beam intensity of two beams which 
IS provides adequate spatial sampling. 

FIGURE 3 illustrates the idealized beam intensity of two more widely 
separated beams which provides spatial samplmg which fails to satisfy the Nyquist 
criterion. 

FIGURE 4 illustrates .the idealized beam intensity of two more widely 
20 separated beams which provides q)atial samplmg which satisfies the Nyquist criterion. 

FIGURE S illustrates the exemplary lobe pattern of an ultrasound beam. 
FIGURE 6 illustrates &e exemplary lobe patterns of two ultrasound 
beams which provide spatial sampling which satisfies the Nyquist criterioiL 

FIGURE 7 illustrates the exemplary lobe patterns of two more widely 
25 separated ultrasound beams which provide spatial sampling which satisfies the Nyquist 
criterion. 

FIGURE 8 illustrates the exemplary lobe patterns of two ultrasound 
beams which provide spatial sampling which fails to satisfies the Nyquist criterion by a 
controlled degree. 

30 FIGURE 9 illustrates an exemplary spatial sampling spectrum. 
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FIGURE 10 illustrates the azimufh and elevation dim^ions of a 
pyramidal volumetric region which is to be eiiicienfiy scanned in accordance wifli the 
principles of file present invention. * 

FIGURE 1 1 illustrates a volumetric ultrasonic diagnostic imaging system 
5 constructed in accordance with the principles of the present mvention. 

FIGURES 12a-12j illustrate the variation in point spread functions at the 
focus of a variety of beams with different combinations of aperture and apodization 
functions. 

FIGURES 13a and 13b illustrate the exemplary lobe pattern of a 
10 relatively narrow ultrasound q)erture in two dimensions with a point spread function 
which is controlled in accordance with the principles of tiie present inventioa 

FIGURES 13c and 13d illustrates the exemplary lobe pattern of a 
relatively broad ultrasound aperture in two dimensions with a pomt spread function 
which is controlled m accordance with the principles of flie present invention. 
15 Referring first to FIGURE 1 an idealized ultrasound beam intensity 

profile 50 is shown. The intensity profile 50 is idealized because it is shown as a square 
function with flie mtensity (amplitude) at a constant maximum intensity and dropping to 
zero intensity on either side of the beam. The abscissa of the beam plot shows that in 
this example the beam extends a half millimeter of distance in azimuth (cross-range 
20 distance) (25.5mm to 26.0mm in fliis example) in the region of the focus of the imaging 
field. 

To adequately spatially sample the imaging field, multiple beams must 
be transmitted which are spaced so as to meet the Nyquist criterion. FIGURE 2 
provides an illustration of a second beam which is transmitted m addition to the beam of 

25 nOURE 1 to adequately spatially sample tiie imagmg field. The second beam has an 
ultrasound beam mtensity profile 52 indicated by the dashed lines. The second beam 
intensity profile is seen to extend from 25.75mm to 26.25mm in this example. Since 
the second beam profile overlaps that of tihe first beam by 50%, the imaging field is 
being spatially sampled so as to meet the Nyquist criterion at this point, which calls for 

30 sampling at twice the frequency of the spatial information. A succession of such beams 
• across the fiiU angular distance of the imagmg field will adequately sample the entire 
imaging field. 
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FIGURE 3 shows the beam intensify profiles 50 and 54 of two beams 
where tfaebeamsaie more widely separated. TTie beam intensity profiles a re ofthe 
same dimensions as in the preceding examples, each extending 0.5mm in azimuth. 
However, in this example the center-to^enter spacing ofthe beams is 1mm in distance 
5 rafter dian the 0.25mm spacing of the preceding example. The wide separation of the 
two beams ftils to satisfy the Nyquist criterion for spatial sampling, and such a beam 
sampling pattem can give rise to the scintillating or shimmering artifects characteristic 
of spatial undersampling. 

In accordance with the principles of the present invention, when the 

10 scannmg beams are more widely separated, the spatial point spread fimctions ofthe 
beams are adjusted to account fi>r the greater center-to-center spacing (reduced output 
line density) of fte beams. As used herein, the point spread function refers to the two- 
way spatial response of a pulse-echo sequence, that is, the beam patterns of a transmit 
beam and its received beam or beams, use for spatial sampling. The point spread 

15 fimction is detennmed by the size of the transducer aperture employed and the 
apodization (weighting or mtensity) function used at the aperture. The drawings herein 
which iUustrate point spread functions generally show a one-way (transmit) relationship 
between the aperture and the pomt spread function at the beam focus. Beam focusmg 
may be layered on top of the aperture control used to define the point spread fimction, 

20 which is generally done by a mechanical lens or electronic delays. FIGURE 4 shows 
two beam intensity profiles 56 and 58 for two beams with a oenter-to-center spacing of 
1mm, the same as the beams of FIGURE 3, but with an aperture fimction that produces 
a broader beam intensity profile (2mm in diis example). It can be seen that the two 
beam intensity profiles 56 and 58 overhp by 50% as in FIGURE 2, resulting m 

25 satisfaction of flie Nyquist criterion for spatial sampling of the imaging area wifli the 
more widely spaced beams. 

The beam intensity profile of an ultrasound beam at the focal plane 
which is transmitted by an array transducer is not square as in the preceding drawings, 
but is more sinusoidal in shape, and due to the finite size ofthe aperture, will generaUy 
30 have a main lobe surrounded by side lobes as shown by the beam intensity profile 60 of 
HGURE 5. Whereas the extent of the beam intensity profiles ofthe precedmg figures is 
clearly delineated by the instantaneous drop to zero at the sides ofthe square profile, an 
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actual beam profile such as the profile 60, which rolls off gradually fix)m its center peak, 
has a spatial extent detennined by the criteria of ttie system designer. One common 
intensity level which is used for the eflBwtive extent of a beam intensity profile is the 
point at which the intensity has rolled oflFby 3 dB firom the intensity peak, indicated by 
5 points 62 and 64 on either side of the main lobe in FIGURE 5. With the 3 dB points 
used m this example, the eiBfective beam dimension for spatial sampling are seen to 
extend over the distance ftom Dl to D2. For adequate Nyquist spatial sampling, the 
3dB point of the adjacent, simUarly dimensioned beam 66 should fell between the 3 dB 
points 62 and 64 of the beam 60, as shown in FIGURE 6. However, if Ae beams are 

10 more widely separated, that is, the widdi of the region being scanned increases or the 
beam density decreases, die point spread functions of die beams are changed so diat the 
3 dB points 72, 78, 74 of the beams 70 and 76 sufficiently over to satisfy the Nyquist 
criterion for spatial sampling as shown in FIGURE 7. 

A point spread function providing a broader main lobe transmit beam ' 

IS will insonify a broader region around the center of the beam profile. This enables tiie 
reception of a gceater number of receive multilines in response to each transmit beam. 
As the transmit beam is broadened, the product of each multiline profile and the 
transmit beam profile provides an improved point spread function for each transmit- 
receive combination. The point spread function in this case is dominated by the 

20 narrower beam profile of each receive multiline. See US Pat 6,494,838 for a 
description of a system which increases flie volumetric line density through multiline 
reception and scanline interpolation. 

Instead of fully satisfying the Nyquist criterion for spatial sampling, it 
may be decided for certain applications to maintain a spatial sampling beam dispersion 

25 which fells short of the Nyquist criterion but nevertheless produces images which are 
satisfactory for tiie given procedure. For example, an obstetrician may be imaging a 
fetus to measure the bones of the fetus for gestational age calculation. In such an exam 
tissue texture may not be important but a higher frame rate may present images of a 
fetus moving m flie womb which may be satisfectorily measured. The obstetrician will 

30 generally be satisfied if the tissue of the anatomical feature is in the correct location, in 
which case a lower spatial ficequency will sufiBce. FIGUllE 8 illustrates two adjacent 
beam profiles 80, 82 which overlap at their adjacent 3dB points 84 (position D2 on the 
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distance axis). While some spatial sampling arti&cts may develop fiom this beam 
spread, they may not be at a level which sigoificanfty impedes Oie ability to make fetal 
bone measurements. If the volume being imaged is increased, tfie aperture of die 
transmit beams may be adjusted to broaden flie beam profiles and hence the extent of 
5 spatial information being interrogated. FIGURE 9 illustrates the relationship between 
the spatial sampling frequency and artifacts developing due to spatial undersampling 
graphically. The area or volume being imaged may be sampled at a spatial samplmg 
frequency f^ which is twice a spatial cutoff frequency fc The anatomical information 
which is being sampled has a band of spatial frequencies 86 which roUs off to an upper 

10 frequency fh. Thus» spatial frequencies above ^ will alias back to a lower frequency fc- 
as indicated by tfie dashed line 88. hi a particular application such aliasing may be 
acceptable; in others^ it may not and spatial sampling fs should be done at a higher 
spatial frequency if textures such as the speckle pattern are desired for the diagnosis. * 
Jn an efiBcient data acquisition design, the sampling bandwiddi or spatial 

IS resolution is matched to &e achievable transducer resolution (which may be 
characterized by the aperture size and die acoustic wavelength) and the desired output 
bandwidth or volume imaging rate. Different combinations of transducer geometry, 
output line density and volume imaging rates lead the efficient design to use variable 
acquisition resolution, hi an ultrasound system with a programmable beamformer, the 

20 spatial pomt spread function can be adjusted to best match the spatial resolution to the 
desired output line density, which will determine the frame rate of the two or three 
dimensional image. In a 3D scanning application where a maxunum volume image rate 
is desired, the point spread function can be altered by adjustment of the apodization of 
the transmit aperture or receive aperture or both to match the sampling resolution to the 

25 line density. A simple example of how this adjustment can be made is iUustrated with 
reference to FIGURE 10. Suppose that the clinician wants to perform 3D imaging of a 
fetal heart Further suppose diat the 3D transducer probe has an array transducer tiiat is 
capable of scanning a pyramidal volume 90 as shown in FIGURE 10. The array 
transducer is located at or just above the apex 92 of the volume 90. Further suppose 

30 fliat the clinician finds that she can capture the entire fetal heart i n a volume which 
measures 30^ in the azimuth direction and 30** in the elevation direction and which 
extends to a depth of 7 cm. as shown in tfie drawing. The round-trip time requured for 
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sound to xeach the 7 cm depth and tetum is assumed to be lOO^sec. in this example. 
This means that the acquisition time for one scanline is lOOfisec. Further assume that 
the clinician desires a fiame rate of 30 volumes per second. From the desired frame rate 
of 30 vol/sec and the line time of lOOjisec/Iine, it is seen that 333 lines can be used to 
5 scan the volume 90 in the time allotted to meet the volume frame rate requirement. 
These lines are to be distributed over the volume 90. Although diflferent line densities 
can be used in the azimuth and elevation directions, in this example it will be assumed 
that a uniform line density in both directions is to be used. The allotted number of Imes 
can be distributed with eighteen lines in the azimuth direction and eigjitem Imes in flie 

10 elevation direction as indicated by the small delineations along tiie base of the volume 
90. For a volumetric sector measuring 30* by 30^ this means that the lines are on 
approximately a 1.6® center-to-center spacing. To meet theNyquist criterion with a 
50% overlap, a pomt spread flmction of 1.6** should be used to satisfy flie Nyquist 
criterion in the elevation and azunutfa directions. In the diagonal durection the volume 

15 will be slightly spatially undersampled, which may be overcome, if desired, by slightly 
widening the beam profile or increasing the line density. The ability to shape tiie point 
spread ftinction in three dimensions with a two dimensional array transducer f urflier 
enables die formation of advantageous shapes of the point spread function. For 
instance, the point spread function can be shaped to yield a hexagonal approximation 

20 for more efficient beam packing in a volume. See, for example, US patents 6,384,516, 
6,497,663, and application serial number 09/908,294, which describe the fabrication 
and use of hexagonal array transducers and beam scanning. 

Thus it is seen diat a method to design the scaiming criteria for a 
volumetric region starts by determining the desned output volume size (30*^ by 30® by 

25 7cm in the above example) and flie desired volume acquisition rate (30 volumes/sec m 
the example). A line density is calculated fliat can be supported by the deshed volume 
size and volume acquisition rate (333 lines/vol. in flie example). The line density may 
be asymmetrical or synwnetrical in all directions. The point spread function is then 
calculated that is required to sample the line density in both azimuth and azimuth (1 .6** 

30 in the example). An apodization function is then chosen that provides the calculated 
point spread function in azunuth and elevation, for the transmit and preferably both the 
transmit and receive beams. An ultrasound system for carrying out this method in 



ATL-335 



accordance wifli the principles of the present invention is shown in FIGURE 11. An 
ultrasonic probe 10 capable of Uiree dimensional imaging includes a two dimensional 
array transducer 12 which transmits beams over a fhree dimensional volume and 
receives single or multiple receive beams in response to each transmit beam. Suitable 
5 two dimensional arrays are described in U.S. patent appl. serial number 09/663,357 and 
in U.S. Patent 6,468,216. The transmit beam characteristics of the array are controlled 
by a beam transmitter 16, which causes the apodized aperture elements of the array to 
emit a focused beam of the desired breadth in a desired direction through a volumetric 
region of fee body. Transmit pulses are coupled from the bdam transmitter 16 to the 

10 elements of fee array by means of a transmit/receive switch 14. The echo signals 
received by fee array elements in response to a transmit beam are coupled to a 
beamfoimer 1 8, where fee echo signals received by fee elements of the array transducer 
are processed to form smgle or multiple receive beams in response to a transmit beam. 
A suitable beamformer for this purpose is described in U.S. patent appl. serial number 

15 09/746,165. Rafeer than housing all of fee bramformer circuitry in fee system 

beamformer 18, fee beamformer circuitry may be distributed between fee probe 10 and 
fee system as described m U.S. patent 6,468,216. 

The receive beams formed by the beamformer 18 are coupled to a signal 
processor which performs functions such as filtering and quadrature demodulation. The 

20 processed receive beams are coupled to a Doppler processor 30 and/or a B mode 
processor 24. The Doppler processor 30 processes fee echo information into Doppler 
power or velocity information. The feree dimensional Doppler information is stored in 
a 3D data memory 32, from which it can be displayed in various formats such as a 3D 
power Doppler display as described m U.S. patmt Re. 36,564. For B mode imagmg fee 

25 receive beams are envelope detected and fee si^ls logarithmically compressed to a 
suitable dynamic range by fee B mode processor 34 and feen stored in fee 3D data 
memory 32. The 3D data memory may comprise any memory device or group of 
memory devices which has three address parameters. The three dimensional image data 
stored in fee 3D data memory 32 nuiy be processed for display in several ways. One 

30 way is to produce multiple 2D planes of fee volume. This is described in U.S. patent 
6,443,896. Such planar images of a volumetric region are produced by a multi-planar 
reformatter 34. The three dimensional image data may also be rendered to form a 3D 
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display by a volume Tenderer 36. Hie resulting images, which may be B mode, Doppler 
or both as described in US patent 5,720.291, are coupled to an image processor 38, fiom 
which they aie displayed on an image displ^ 40. 

Tn accordance wifli flie principles of flie present invention the ultrasound 
5 system of FIGURE 1 1 includes a beamformer controller 22 which controls both the 
beam transmitter 16 and die receive beamformer 1 8. The beamformer controller 22 is 
responsive to a user interface 20 by which a clinician may set imagmg parameters for 
the beamformer confroUer. The clinician may input values for the azimuth and 
elevation widths of a volumetric scan region, the depth of the scan region, and a 

10 required fiame rate, for instance. Ultrasound systems such as ttiose available fiom 
Philips Ultrasound lie, can choose fliese initial parameter settings automatically m 
response to the selection of an exam type by die clinician, a feature known as 'Tissue 
Specific Imaging." From these parameters tiie beamfomer controller can calculate tite 
number of lines which can be used to scan the votometric region and the line density as 

15 discussed above, and the point spread function needed for that line density. Since tiie 
focal plane point spread function is the Fourier transform of the aperture function, die 
beamformer controller 22 can perform an inverse Fourier transform of the point spread 
function to calculate die needed anay aperture. Alternately, the parameters for the 
desired point spread function can be precalculated and stored on the system for 

20 implementation togeflierwifli die programmed focus parameters. Itmayalrobe 

sufBcient to determine the point spread fimction "on die fly" by dioosing an appropriate 
aperture, as the point spread fimction is approximately inversely proportional to die 
aperture fimction. As die signals to or fiom die transducer elements of flie aperture arc 
shaded (differently weighted, or apodized), die point spread fimction will broaden to 

25 accommodate a greater line spacing (lesser line density). Stated anoflier way, die beam 
widfli is inversely proportional to die aperture widdi. By varying die number of 
transducer elements and dieir locations of die active aperture for transmit and/or 
receive, and die weightings of die signals to or from fliose elements (which also affects 
side lobe characteristics), die widtii of die mam lobe of die acoustic beam is tailored for 

30 die desired point spread function. See Optics. Second Editinn hy Bng^nft Ut^ht 
(Addison-Wesley Pub. Co.) at C3i. 11, and Introduction To Fourier Op rin!^ by J. W. 
Goodman (McGraw-HiU Book Co.) at Ch. 4, in which diese princq)Ies are illustrated in 



ATW35 



the field of optics. 

FIGURES 12a-12j illustrate the variation in point spread functions with 
different aperture and apodization combinations for volumetric imaging in accordance 
with the principles of the present invention. In each of these drawings the numbers at 
5 the base grid refer to size measures in the elevation and azimuth directions. For an 
anay of transducer elements that are uniformly sized and spaced m both elevation and 
azimuth, the base grid of these drawings would correspond to the elements of a 64 
element by 64 element transducer array. The height of the beam pattern above each 
point on the grid (element) conresponds to the relative apodization function at tfiat 

10 particular point (element of flie array.) Thus, the shape of grid area beneatii each beam 
pattern mdicates the elements used for the active aperture and tiie shape of the beam 
pattern above fliose elements shows the apodization function used to produce the point 
spread function at the focus. In FIGURE 12a the active aperture comprises a 
symmetrical central area of sixteen elements in azimuth and sixteen elements in 

15 elevation. A Banning window is used for apodization in bofli flie elevation and azimuth 
directions as shown by the shape 100. This aperhire function will produce a point 
spread function or beam pattern 102 at the focus as shown in FIGURE 12b wifli the 
greatest intensity (greatest weighting) in the center and declining smoothly and 
uniformly in both the elevation and azimufli directions from the center. The Harming 

20 window apodization results in relatively low side lobe levels. 

FIGURE 12c illustrates an ap^ture function 1 10 produced by an 
asymmetrical 1:2 aperture of sixteen elements in azhnufii and fliirty-two elements m 
elevation. A Hanning window is used to smoothly apodize the aperture in each 
dimension from a common central point in flie center of the transducer. This q)erture 

25 function produces fihe point spread function or beam pattern 1 12 as shown in FIGURE 
12d. The aperture function which is broader in tfie elevation dhnension is seen to 
produce a point spread fimction 1 12 which is narrower in the elevation dimension at the 
focus. A point spread functions such as fliat shown in FIGURE 12d might be used 
when a greater spatial resolution or different number of multilines is desfred in one 

30 dimension as opposed to die other. 

FIGURE 12e shows a reversal of the aperture function of FIGURE 12c. 
In this case the aperture function 120 has a greater breadth in the azimuth dimension. 
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producing a beam pattern or point spread fimction 122 which is narrower in flie aamuth 
dimension, as shown in FIGURE 12f. This pomt spread fimction might be used when 
greater lateral resolution in the azimuth dimension or higher multiline order in the 
elevation dimension are desired. 
S FIGURE 12g illustrates fhe aperture function 130 of a 1 :2 aperture with 

unvaiying (rectangular) apodization. The lack of a smootii apodization function 
produces a beam pattern or point spread function at the focus which exhibite a main 
lobe 132 and side lobes 134 m both the elevation and azimuth dimensions. If a 
smoo^y varymg Harming window were used for the apodization function in tiie 
10 elevation dimension as shown by the apertuie function 140 in FIGURE 12i, tiie 
resultant point spread function 142 would have tiie significant side lobes 144 in tiie 
azimuth dimension but not the elevation dimension as shown in FIGURE 12j. 

FIGURES 1 3a-l 3d illustrate how the aperture function can be changed 
by the setting of the s^erture and apodization functions by the beamformer controller to 
1 5 produce a broader or narrower pomt spread fimction that provides the desired spatial 
samplii^ fiequency. FIGURE 13a shows an asymmetric three dimensional aperture 
fimction 150 witii an active aperture of eight elements by sixteen elements and Harnimg 
window apodization in both elevation and azunutii. This aperture fimction produces -a 
point spread fimction 1 52 at tiie focus as shown in FIGURE 13b. The point spread 
20 function 152 is relatively narrow m tiie elevation dimension and broader in the azimutii 
dimension with relatively low side lobe levels. If the volume which is scanned with 
beams of Ibis nature is to be scanned at a higher fi:ame rate, an aperture function 160 as 
shown in FIGURE 13c could be used. As FIGURE 13c shows, the new aperture 
function only occupies an aperture of five elements by eight elements and is apodized 
25 ' wifli a Hanning window. This aperture fimction will produce a much broader point 
spread function 162 at the focus as shown by FIGURE 13d. It can be seen that fewer 
beams of tiie beam pattern of FIGURE 13d will be needed to scan a volume of a &vm 
size tiian beams of the beam pattern of FIGURE \3h, thus enabling flie volume to be 
scanned at a higher volume display rate. 
30 An embodiment of tiie preset invention can, if desu:ed, advantageously 

provide increased scan deptiis as &e point spread function is varied. The acoustic 
output of medical ultrasound transducers is regulated in most countries by maximum 
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aUowable levels of peak acoustic pressure and of average or long-tenn thennal energy, 
to the United States fliese parameters are controlled by limiting the Mechanical iidex 
and IsPTA of the acoustic transmissions. FIGURE 13b illustrates the beam profile of a 
relatively narrow point spread function where most of the energy of the transmitted 
5 beam is concentrated in a relatively narrow central lobe 152 which extends over a 
relatively nanow central area of the array and hence is relatively concentrated in tiie 
body. To avoid exceeding peak acoustic pressure Ihnits the energy in the relatively 
tightiy contained area of the central lobe 152 must be limited to relatively low levels 
and the narrow lateral extent of the beam profile limits the overall energy provided by 
10 thebeam. FIGURE 13d, on the other hand, iUustrates the beam profile of a relatively 
broader point spread function which may be used when the clinician calls for a higher 
volumetric fiame rate or a wider volumetric region. In such cases broader pomt spread 
fimctions for a decreased beam density meeting the Nyquist or Nyquist-related criterion 
are employed. For tiiis beam the energy fiom the array transducer is distributed over a 
15 greaterareainthebody,theareaoffliebroaderbeampatteml62. More energy can be 
transmitted by fewer transducers smce fee point spread function exhibits this broader 
lobe. Consequently the transmitted beam contains more energy and can penetrate to 
greater depths in the body, returning clinically useful echo information from greater 
deptiis without violating acoustic output limits. Accordingly, by varymg die total 
20 acoustic ou^ut power in concert with changes in the point spread fimction. the changes 
in the point spread function can advantageously be used to increase acoustic penetration 
and die clinically useful depth of die image. 

As the point spread function is relaxed (broadened), the effective focal 
range of the beam extends over a wider range of deptiis. The extended depth of focus 
25 meansthatanincreaseddepthoffieldcanbeimagedandremainmfocus. Anincreased 
depfli of field can reduce the need for multiple focal zones, fliereby increasing die 
vohmietric fiame rate. A reduction m die need for multiple focal zones is very 
significant in fliree dimensional imaging because tiie volume fiame rate reduction 
caused by multiple transmit focal zones can be severe. 

considerations may also aflfect the design of the apodization 
fimction. For instance, a phased array which is angularly steered will perform 
differenUy at die sides of the array where die steeply steered beams cause transducer 
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acceptance angle effects. When the angular samplmg densi^ is to be maintained 
constant throughout the volume, the apodization function may vary with beam angle to 
compensate for transducer acceptance angle effects that otfieiwise would lead to a 
variable point spread function in different parts of die image region. 

5 Other variations will readily occur to those skilled in the art For 

instance, the ability to shape the point spread function enables the beam density and 
beamwidth to be varied throughout the image field. A higher beam density could be 
employed in the center of a volume, with a relaxed point spread function and lower 
beam density used at tiie lateral extremes of the volume. The beam density can be 

10 varied continuously from the center to the sides of the volume being scanned. 

An embodiment of tfie pr^ent invention can be used as desired to 
improve the mformation content of the echo mformation and the information movement 
efiBciency by not acquiring more resolution than can be utilized. It can also provide a 
more optimal sampling fimction by usmg die ^erture function to Kmit die spatial 

15 (azimudi and elevation) bandwiddi for duee dimenaonal imagmg. 
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WHAT IS CLAIMED IS: 

1 . An ultrasonic diagnostic imaging system for duee dimensional 
scanning conqnising: 

5 an array transducer having a plurality of transducer elements; 

a beamformer coupled to the array transducer which causes flie 
transducer to scan a volumetric region with a plurality of hansmit beams and to receive 
echo information in response to transmit beams, the beamformer controlling the point 
spread functions of beams transmitted and/or received by Ae beamfoimer, 

M> iniage processor coupled to the beamformer which produces image 
signals in response to the echo information; and 

a display coupled to fbe image processor, 

wherein beams produced by the beamformer exhibit a first point spread 
fimction when the volumetric region is scanned with a first line density and a second 
15 point spread fimction when «ie volumetric region is scanned with a second line density. 

2. The ultrasonic diagnostic imaging system of Claim 1, wherein the 
point spread fimction comprises the two-way spatial response at a focal region of pulse- 
echo spatial sampling of ihe volumetric region. 

20 

3. The ultrasonic diagnostic imaging system of Claim 1, wherein die 
transmit beams exhibit a relatively narrower beam profile at the focus when scanning 
the volumetric region with a first line density, and the transmit beams exhibit a 
relatively broader beam profile at the focus when scanning the volumetric region with a 

25 second line density which is less flianfte first line density. 

4. The ultrasonic diagnostic imaging system of Claim 3, wherein 
adjacent beams overlkp at substantially the same intensity levels when scanning the 
volumetric region witti the first and second line densities. 

30 

5. The ultrasonic diagnostic imaging system of Claim 4, wherein fee 
transmit beams satisfy fee Nyquist criterion for spatial sampling of fee volumetric 
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region to substantially the same degree. 

6. The ultrasonic diagnostic imaging system of Claim 1, wherein the 
point spread functions satisfy the Nyquist criterion for spatial sampling of flie 
volumetric region to substantially the same degree. 

7. The ultrasonic diagnostic miaging system of Claim 1, wherein flie 
beam point spread function exhibit both an azimuth dimension and an elevation 
dimension; 

wherein pomt spread function is symmetrical in both flie azimuth and 
elevation dimensions. 

8. The ultrasonic diagnostic imaging system of Claim 1, wherein flie 
beam point spread function exhibit bofli an azimufli dimension and an elevation 
dimension; 

wherein point spread function is asymmetrical in the azimuth and a 
elevation dimensions. 

9. An ultrasonic diagnostic imaging system for three dimensional 
scanning comprising: 

an array transducer having a plurality of transducer elements; 

a beamfonner coupled to the array transducer which causes the 
transducer to scan a volumetric region wifli a plurality of transmit beams and to receive 
echo information in response to transmit beams, the beamfonner controlling the point 
spread functions of beams transmitted and/or received by the beamformer by control of 
the aperture function of flie array transduce^ 

an image processor coupled to flie beamfonner which produces image 
signals in response to flie echo information; and 

a display coupled to the image processor^ 

wherein flie beamformer utilizes a first aperture function when the 
volumetric region is scaimed with a first line density and a second aperture fimction 
when the volumetric region is scanned wifli a second line density. 
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10. The ultrasonic diagnostic imaging syst^ of Claim 9, wherein flie 
aperture fonction comprises the combination of the elemmts used in an active aperture 
of the anay transducer and the apodization function of the elements of the active 
aperture. 

11. The ultrasonic diagnostic imaging system of Claim 10, wherein 
the apodization function is controlled to match the point spread function to the line 
spacing when scanning the volumetric region with the first and second line densities. 

12. The ultrasonic diagnostic unaging system of Claim 1 1, wherein 
the first line density is greater dian the second Ime density; and 

wherein the apodization function is controlled to scan an increased 
depth-of-field when scanning Ihe volumetric region with the second line density. 

13. The ultrasonic diagnostic imaging system of Claim 10, wherein 
the apodization function comprises the relative weighting of signals of the respective 
elements of flie active aperture durmg a transmission or reception event. 

14. The ultrasonic diagnostic unaging system of Claim 9, wherem the 
first and second aperture functions satisfy the Nyquist criterion for spatial sampling of 
the volumetric region to substantially the same degiee. 

15. The ultrasonic diagnostic imaging system of Claim 14, herein 
the first and second aperture functions both substantially exactly satisfy the Nyquist 
criterion for spatial samplmg of flie volumetric region. 

16. The ultrasonic diagnostic imaging system of Claim 10, wherein 
the scanning beams exhibit a substantially constant angular sampling density; and 

wherein the apodization function is varied as a function of beam angle to 
compensate for transducer acceptance angle effects. 
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17. Inm ultrasonic diagnostic imaging system for volumetric 
scanning and which includes a user interface, a mettiod for detennining the point spiead 
function used to spatially sample a volumetric region comprising: 

determining the desired size of tfie volumetric region to be scanned; 
5 determining the desired volume acquisition rate; 

calculating the line density for scanning the volumetric region of the 
desired size at the desired volume acquisition rate; and 

calculating the point spread function which will spatially sample the 
volumetric region at the line density. 

10 

18. The method of Claim 17» wherein calculating die point spread 
function further comprises calculating the point spread function which satisfies the 
Nyquist criterion for spatial sampling of the volumetric region to a desired degree. 

15 19. The method of Chiim 17, further comprising: 

determining an aperture function tiiat provides the calculated point 
spread function. 

20. The method of Claim 1 9, wherein determining an aperture 
20 fimction comprises determining an apodization function for an active aperture that 

provides tiie calculated point spread function. 

2 1 . The method of Claim 1 7, wherein detennining the desired 
volume acquisition rate comprises determining the volume fi:ame rate of display. 

25 

22. The method of Claim 1 7, wherein calculating the point spread 
function comprises detennining an aperture function which is approximately inversely 
proportional to a desired point spread function. 
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IMPROVED ULTRASONIC VOLUMETRIC IMAGING BY 
COORDINATION OF ACOUSTIC SAMPLING RESOLUTION, 
VOLUMETRIC LINE DENSITY, AND VOLUME IMAGING RATE 

5 

Abstract of the disclosure: 

In an ultrasonic diagnostic imaging system which scans a volumetric 
region, the sampling bandwidth or spatial resolution is matched to the achievable 

10 transducer resolution detennined by the aperture size and the wavelength and the 

desired output bandwidth or volume imaging rate. In an illustrated embodiment this is 
done by controlling fiie spatial point spread function of ttie beams used to scan the 
volumetric region to provide a mote optimal relationship between the acoustic sampling 
resolution, the desired ou^ut line density, and the volume imaging rate. The benefits of 

1 5 this optimization can be to maximize the information content and information 
movement efficiency by not acquiring more resolution tiian can be utilized, and to 
provide a more optimal sampling fimction, using the aperture fiinction to limit the 
spatial bandwidth. 
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